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Abstract. In this second paper devoted to the equatorial symmetry/antisymmetry ol 
■ stationary axisymmetric electrovac spacetimes we show how two theorems proved in our 

previous paper (Ernst, Manko and Ruiz 2006 Class. Quantum Grav. 23 4945) can be 

(N ■ 

utilized to construct exact solutions that are equatorially symmetric or antisymmetric. 
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1. Introduction 



Theorems which allow one to single out "equatorially symmetric" or "equatorially 
antisymmetric" spacetimes from the totality of stationary axisymmetric electrovac 
solutions of the Einstein-Maxwell equations were recently published p] by Ernst, Manko 



and Ruiz. Specifically, in terms of the complex Ernst potentials [2] £ and $, equatorially 
symmetric spacetimes satisfy the conditions 

£(-z,p)=£(z,p)*, $(-z,p) = e 2i5 $(z,p)*, 6 = const, (1) 

while equatorially antisymmetric spacetimes satisfy the conditions 

£{-z, p) = £(z, p), $(-z, p) = ±$(z, p). (2) 

Alternatively, in terms of the axis data e + (z) and f+(z) that comprise the values of the 
potentials £ and $ on the upper part of the symmetry axis, analytically extended 
to negative values of z, asymptotically flat stationary axisymmetric solutions with 
equatorial symmetry were shown to satisfy the relations 

e + (z)[e + (-z)Y = l, f + (z) = -e 2iS if + (-z)Ye + (z), (3) 

while asymptotically NUT equatorially antisymmetric solutions were shown to satisfy 
the relations 

e+(z)e+(-z) = 1, f + (z) = Tf+(-z)e+(z). (4) 
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A large number of astrophysically relevant electrovac solutions were considered in 
the papers [3], H] and [5j [6]; these correspond to rational axis data, which for some 
sufficiently large value of N can be expressed in the form 

z N + Ek=ihz N ~ k R{z) 
ELiC k z N - k Q(z) 
z N + ELihz N ~ k R{z) 
where a k , b k , c k are arbitrary complex parameters. By systematically exploiting the 
definition ([I]) and the condition (j3J), we shall in this paper construct new equatorially 
symmetric solutions corresponding to such rational axis data. Similarly, using ([2]) 
and (HI), we shall construct new equatorially antisymmetric solutions. Finally, we 
shall discuss the prospects for constructing all equatorially symmetric or antisymmetric 
solutions that correspond to data (15al) and (15bjl with arbitrarily large values of N. In 
particular, we shall show how this can be accomplished for vacuum fields. 

2. Axis data of equatorially symmetric/antisymmetric fields 

To single out equatorially symmetric axis data, we first substitute e + (z) from ([5a]) into 
the first condition from (J3j), whence we obtain 

P(z) = (-l) N R*(-z). (6) 

From this we obtain the first restriction upon the parameters; namely, 

a k = (-!)%, k = l,...,N. (7) 

The substitution of e + (z) and f+(z) from (!5a|) and (|5b| into the second condition from 
(J3]) then leads, after taking account of (jSj), to 

Q( z )+e 2i5 (-l) N Q*(-z)=0, (8) 

which imposes a second restriction upon the parameters; namely, 

c k + e 2i5 (-l) k cl = 0, k = l,...,N. (9) 



Turning now to the equatorially antisymmetric case, we first get from (l5al) . ( I5bl) 
and (jlj) the relations 

P(z) = (-l) N R(-z) and Q(z) =F (-l) N+1 Q(-z) = 0, (10) 

and then obtain from ( TlOl the following restrictions upon the parameters: 

a k = (-l) k b k , c k T(-l) k+1 c k = 0, k = l,...,N. (11) 

Formulae ([7j), © and ffTTj) . which provide the general forms of the axis data fl5aj) 
and ( I5bl) in the equatorially symmetric/antisymmetric cases, will simplify the search 
for examples of axis data from which may be constructed Ernst potentials £ and $ of 
spacetimes of particular physical interest. 
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As an illustration of the use of the above relations let us take N = 2 in the axis data 
f l5a|) and f!5bj) . From ([7]) we find that in the equatorially symmetric case the function 
e + (z) has the form 



e+{z) 



b\z + b* 2 



(12) 



z 2 + b x z + b 2 
while from ([9]) we see that f+(z) has the form 

e iS (qz + ic) 

J + \ z ) ~ i , u nr> 
z 2 + biz + b 2 

where q and c are arbitrary real constants. On the other hand, in the equatorially 
antisymmetric case, from ffTTj) one finds that 

z 2 - biz + b 2 



(13) 



while for / 4 



and 



(2;) one arrives at the following two possibilities: 

c 2 



.A 



2 2 + b\z + 6 2 



(14) 



(15) 



(16) 



z 2 + b 1 z + b 2 

Note that all parameters in formulae ()14p - (jl6p may be assigned arbitrary complex 
values. 

In the next section we shall determine the Ernst potentials S(z,p) and 
corresponding to three specifications of the axis data; namely, that specified by ([121) 
and ( fTBl . that specified by ( TT4"1) and (EES]), and that specified by ( 1141 and (|T6l) . 



3. £ and $ for rational axis data 

The general solution arising from the rational axis data (I5ap and (15bl) is known [5], 
but the description of the equatorially symmetric and antisymmetric subfamilies of 
this solution in terms of the parameters it contains has not yet been treated in any 
publication. Recall that 



where 



£ = E + /E„ 
1 

±1 
±1 

£4 



$ = F/E. 



1 


1 


n 


r 2 N 




&2N - Pi 


r% 


r 2N 


oil - Pn 


ot 2 N — Pn 


hx(ai) 
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(17c) 



hi{a n ) 



ei 



p 2 + (z- a n ) 2 , 

e r + 2/;/K), 

2n^i(A-« r 



A? 



/K) = E 



./'/ 



a. 



i=i 
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- Pi 
fifk 



and 

' (17d) 

(17e) 

n^(A - A) nl v =1 (A - " fe ti A - pi ' (17f) 

Here E± and F are determinants of order 2N + 1, the arbitrary parameters are Pi, a n , 
fl, and all other constant objects are expressible in terms of these parameters. 

The constants Pi can take on arbitrary complex values, while the a n take on 
arbitrary real values or occur in complex conjugate pairs. The parameters ft are 
arbitrary complex constants, and together with the Pi they enter the axis data 



N 



e + (z) 



i + E 



ei 



N 



ft 



i=i - Pi i=i z ~ Pi 

which is a different but mathematically equivalent representation of the data (15al) and 
f l5b|) . The two sets of the parameters, {Pi,oi n ,ft} and {Pi,ei,ft}, are mathematically 
equivalent too, and they are related to each other via the algebraic equation 

S(z) = e + (z) + e\{z) + 2f + (z)f* + (z) = 0, (19) 

which is satisfied identically by the parameters a n ; i.e., S(a n ) = 0, n = 1, 



,2JV. 



3.1. Axis data l[TE) and H 



The axis data (fT2ll and (fl~3l) involves (in addition to the real constant 5 which we shall 
suppress) two complex parameters {61,62} and two real parameters {q,c}. We shall 
express the two complex parameters in terms of four real parameters {m, a, 6, k} that 
are of more direct physical significance. The axis data then assume the form 

ia)(z + ib) — k 



m 



e+{z) 



/+(*) 



[z + m — ia) [z + i6) — k ' 
qz + ic 
(z + m — ia) (z + ib) — k 



(20a) 
(20b) 
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The solution corresponding to such axis data should be useful for modelling the exterior 
field of a neutron star, for the mass-quadrupole and angular-moment um-oct up ole 
moments it possesses are arbitrary and independent. This is expected to result in a 
generalization of two exact solutions [H [8] which were recently shown to be useful in 
this way [HI [TQl [TT]. Incidentally, there were various typographical errors in [T2] that 
will be corrected in the present paper. 

The first step in obtaining the potentials S and $ from the axis values (I20aj) and 
( I20bl) consists in getting the corresponding expressions for a n in terms of the parameters 
m, a, b, k, q, c. The substitution of (I20ap and (120bl) into (|T9|) yields the biquadratic 
equation 

z 4 - (m 2 - a 2 - b 2 - q 2 + 2k)z 2 + (k — ab) 2 - b 2 m 2 + c 2 = 0, (21) 
whence we get 

Oil = -«4 = ^(k+ + K-), "2 = -«3 = \{ K + - «-), (22) 

where 

re± = y/m 2 - a 2 - 6 2 - g 2 + 2A; ± 2d, d = ^/(fc - a&) 2 - b 2 m 2 + c 2 , (23) 
thus also identifying four functions r n : 

ri = -R_, r 2 = r_, r 3 = r + , r 4 = R + , (24) 

where 



= vV + ± U K + + K -)] 2 ' r ± = \A° 2 + [* ± K K + - ( 25 ) 
The next step is obtaining the expressions for e« and //. It is easy to see that 
the axis data (120aj) and f!20bj) can be rewritten in the form 

Z — Pl 2 — p 2 2 — Pi ^ - P2 

where 

/3 1 = -i[m-i(a-6) + v / ^D], & = ~[m - i(a - 6) - «/D], (26b) 
2m(A + ib) 2m{(3 2 + ib) , 

ei = ^^T' e2 = ^^r' (26c) 

/: f^, / 2 = f^ ^=[— Ka + & )] 2 + 4., (26d) 

Pi — P2 P2 — Pi 

so that now all the quantities entering the determinants (117bj) and fllTcj) are defined in 
terms of the parameters of the axis data (I20al) and ( I20bl) when we take into account 
that the f(a n ) are obtainable from f+(z) by simply changing z to a n , 

/(an) = 7 : -r-, ■ T7T r, (27) 

(a n + m — ia)(a n + ib) — k 

while hi(a n ) are combinations of ei, fi and f(a n ). 

The final step is the expansion of the determinants (117bl) and (I17cl) . which in the 
case N = 2 have the order 5, followed by the substitution of e/, //, f(a n ) and 
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hi(a n ) as given by the formulae fT2"2"j) . (1211) . (I26a,j) (126d[) and (12"T|) . The evaluation and 
simplification of the resulting expressions can be carried out with the aid of a symbolic 
computer programme. The final formulae for 8 and $ ar^jj] 

£ = 4— $ = ^— , (28a) 



where 



A = d[m 2 -q 2 - (a - fo) 2 ][/4(i2 + r- + + k 2 _ (R + r + + i2_r_)] 

+ [(a - fo) 2 (afo - k) + (a - 6) (eg - 6m 2 ) + m 2 fo 2 - c 2 ] 

x [K5.(i2+r_ + i?_r+) - K 2 _(R + r + + i?_r_)] 

+ 4d[k(m 2 - g 2 ) + (c + aq)(bq — c)](i2+i2_ + r + r_) 

+ iK + K_{/s; + [(a — fo)(afo — k — d) — m 2 b + cg](i? + r_ — R-r + ) 

+ /c_[(a - 6)(a6 - A; + d) - m 2 b + cg](i?_r_ - i?+r + )}, (28b) 
Z? = mK + K-{d[K+K-(R + + i?_ + r + + r_) 

- (m 2 - a 2 + fo 2 - g 2 )(i?+ + iL. - r + - r_)] 

+ i6d[(«+ + «;_)(#+ - R-) + (k+ - K_)(r_ - r+)] 

+ i(m 2 b — cq + ak — a 2 b) 

x [(«+ + K _)(r + - r_) + (k+ - K_)(i2_ - #+)]}, (28c) 
C = K+K _{d[g(m 2 - (a - fo) 2 - g 2 ) - 2c(a - fo)](r + + r_ - i? + - i?_) 
+ dqn + K_(R + + R_+r + + r_) + icd[(K + + «_)(.&+. - 
+ (« + — «_)(r_ — r + )] + i[(fc — afo)(c + ag — fog) + q{m 2 b — eg)] 
x [(«+ + «;_)(r + - r_) + («+ - /c_)(i2_ - R+)}}. (28d) 



3.2. ylxis data (fl^p and < f73]) 

Proceeding in a similar way, the axis data (I14R and fllSH can be reexpressed in the form 



/N z - k - m - i(a + v) z + k — m + i(a — v) 

e +\ z ) = i ~, r ■ ; ~, r, (29a) 

z — k + m — i(a — z/j z + k + m + i(a + i/j 

/+(*) = 7-^ v % t^, v vP ( 29b ) 

where m, a, u, k, q and b are arbitrary real parameters. Expressions (I29ap and ( I29bl ) 
generalize the axis data of the Breton-Manko solution [13] and can be interpreted as 
representing two identical counter-rotating electrically and magnetically charged masses 
possessing a NUT parameter, m, q and fo being, respectively, the individual mass, electric 
and magnetic charges, a the rotational parameter, v the NUT parameter and k the 
separation constant. 

|| The expression for A in p~2] (see p. 3068) was published with some misprints. However, all the rest 
of the formulae determining the metric functions (3.11) of [12] are free of misprints, albeit there is a 
missing right square bracket in the second line of the function D just after the first curly bracket. 
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Repeating now the steps outlined in the first example, this time starting with 
the axis data (129aj) and f!29bj) . one finally arrives at the following expressions for the 
potentials £ and $: 

A + B A + B v ' 



where 



and 



A = (m 2 + v 2 - q 2 - b 2 ){[(m 2 + k 2 + v 2 + a 2 ) 2 - A(mk + av) 2 } 
x (R + - i?_)(r+ - r_) - 5a + a^(R + + R-){r + + r_)} 
+ 2« + a_{[(m 2 + v 2 - q 2 - b 2 ){m 2 - k 2 + a 2 - v 2 ) + 2(mu - ka) 2 } 
x (R + R_ + r + r_) + 2id(mz/ - ka)(R + R^ - r + r_)}, (30b) 

i? = Ada + a-{m + iz/){[m 2 + v 2 — q 2 — b 2 + \{mv — ka)} 

x (R + + R_ - r+ - r_) - d(i?+ + R- + r + + r_)}, (30c) 

C = 4da + a_(q + ib){[m 2 + is 2 - q 2 - b 2 + \{mv - ka)} 
x (R + + R_ - r + - r_) - + i?_ + r + + r_)} 
= (g + i6)S/(m + iz/), (30d) 

R± = \jp 2 + (z ± ct + ) 2 , r± = y^ 2 + (2 ± a_) 2 , (30e) 
a± = V5 ± 2d, 5 = m 2 + £; 2 + 3z/ 2 -a 2 -2(g 2 + & 2 ), (30f) 



d = y (m 2 + z/ 2 — g 2 — b 2 )(k 2 + 2z/ 2 — a 2 — g 2 — 6 2 ) — (my — A;a) 2 . (30g) 
It is easy to verify that the Ernst potentials given by (1 3 ah ( [30g] ) satisfy the 



conditions 

£(-z, p) = £(z, p), $(-z, p) = $(z, p). (31) 

It should also be pointed out that the solution so defined is not asymptotically flat due 
to the presence of the non-zero NUT parameter v. When v — b — 0, it reduces to the 
Breton-Manko electrovac solution [13J. 

3.3. Axis data (fT^j and (lb]) 

We can reexpress the axis data ( |T4l) and ( {TBI in the form 

z — k — m — i(a + v) z + k — m + i(a — v) 

e +( z > = — r~i ~i v ' i / i — ^7— — 7' ( 32a ) 

z — k + m — i{a — v) z + k + m + i[a + v) 

, , s = 2(x + ic) 2b 

^' [z-k + m-\{a-v)][z + k + m + \{a + v)Y ^ ' 

where the four parameters m, a, v and k have the same meaning as in the previous 
example, while \ an d c describe, respectively, the electric-dipole and magnetic-dipole 
moments of the sources. 
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The corresponding expressions for the potentials S and $, obtainable from the 
determinants (I17bl) and (|17c| . have the form 

A-B C , . 

£ = -A n> $ = "1 n> 33a 

A + S A + S V 7 



where 



A = {5[( m 2 + z/ 2 ) 2 + ( mi/ - kaf] - d 2 (3m 2 -k 2 + u 2 + a 2 )} 
x (R + - i?_)(r+ - r_) - a+a_[5(m 2 + z/ 2 ) - % 2 - c 2 ] 
x + r_)(r + + r_) + 2a + a_{[(m 2 + v 2 ) 2 + (mi/ - ka) 2 - d 2 } 
x (R + R_ + r + r_) + 2id(mz/ - ka){R + R„ - r + r_)}, (33b) 

B = 4da + a_(m + iz/){[m 2 + v 2 + i(mz/ — fca)] 

x (i?+ + - r + - r_) - d(R + + £L + r + + r_)}, (33c) 

C = 4d(x + ic){[m 2 + z/ 2 - i(mi/ - fca)] - i? + ) - a + (r_ - r+)] 

+ d[a_(#_ -R+) + a+(r_ -r+)]}, (33d) 



and 



R± = VP 2 + (z ± «+) 2 , r± = VP 2 + (* ± «-) 2 , (33e) 



a± = V6 ± 2d, 5 = m 2 + k 2 + 3v 2 - a 2 , (33f) 



d = y ("^ 2 + v 2 ){k 2 + 2z/ 2 — a 2 ) — (mi/ — A;a) 2 — x 2 — c 2 . (33g) 

Please note that the 5 defined in ( 133f|) should not be confused with the phase 5 employed 
in the symmetric case! 

The Ernst potentials given by (j33aj) -( |3"3"g"| ) satisfy the conditions 

£{-z, p) = £{z, p), p) = -$(*, p), (34) 

and the solution describes two counter-rotating masses endowed with the NUT 
parameter and with both electric and magnetic dipole moments. When v = 0, the 
solution is asymptotically flat. 

4. Prospects for a completely general solution 

A similar procedure to that which we employed for N = 2 should work also for N = 3 
and N = 4, but the question arises whether or not one can solve the general equatorially 
symmetric/antisymmetric problem for arbitrarily large values of N. At the present time, 
we do not have an answer to this question, but in this section we shall show that in the 
vacuum case this is indeed possible, which gives us hope that the same will ultimately 
be possible in the case of non-vacuum electrovac fields. 
Instead of the parameters the quantities 

X n ss X(a n ) = ffi"""^ , (35) 

11^=1 (On - Pi) 
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can be used in the vacuum solution, in which case the formulae for 6 assume the form 
P A + r (oa \ 



where 



A= J2 K-iN r h- ■ - r i N , (36b) 

ii<...<ijv 

r= X] 7ii...ijv-i r «i • ' • T iN-li (36c) 



«1<...<IJV_1 



and 



N 

\ ■ 



K...i N = {-l) n+ - +tN V{a ill ... 1 a iN )V{a ilil ... 1 a t , N )X[X l . 

n=l 

{i' x < ... <i r N ,i'^i) (36d) 

N-l 

7 n .., N _ 1 = (-ir + --- +l ^V(a ll ,...,a iN _jV(a^ J] 

n=l 

(i[ < . . . < i' N+1 ,i' ^ i) (36e) 
V(a 1 ,...,a k ) = n (ai-aj). (36f) 

l<i<j<k 

In order to show how the subclasses of the equatorially symmetric or antisymmetric 
solutions are contained in (I36al) - (136fl) we must find 2N restrictions on the parameters 
a n and X n which would be equivalent to iV conditions ([7]) or ffTTj) on the complex 
constants a k - The conditions on the parameters a n can be trivially established since 
these parameters determine the location of sources on the symmetry axis. Assuming 
without loss of generality the order 

Re («i) > Re (a 2 ) > . . . > Re («2at-i) > Re («2jv) (37) 

we immediately arrive at N relations between the parameters a n . These relations are 
the same in the equatorially symmetric and antisymmetric cases; namely, 

«1 = -OC2N, «2 = -OC2N-1, ■ ■ ■ , «7V = -«JV+1- (38) 

A simpler way to write this is 

afc = -a 2N -k+i, k — 1, . . . ,N. (39) 

The above conditions determine the equatorially symmetric/antisymmetric character of 
the mass distribution. 

The restrictions on the quantities X n can be obtained from the vacuum 
specialization of equation (19); namely 

e+(z) + e* + (z) = (40) 

with the aid of the conditions (Q and ffTUl) . In the equatorially symmetric case the 
substitution of 

f~l z - Pi R{z) 
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into (HP]) gives, after taking into account ([H]), 
P(z) P*(z) (-1) N 

Qrx + —t^t = / \ / J R*(z)R*(-z) + i?(^)-R(-^r)] = 0, (42) 
#*(z) R(z)R*(z) [ v 7 v ; v ; v ;j , v ; 

whence, bearing in mind that a n satisfy the above equation identically, we get 

R*(a n )R*(-a n ) , . 

— — — + 1 = X(a n )X{-a n ) = -1. (43) 

R{a n )R{-a n ) 

Therefore, taking into account (1391) . we arrive at the remaining N restrictions on the 
parameters X n ; namely, 

X k X 2N _ k+1 = -1, k = l,...,N. (44) 

In the equatorially antisymmetric case we proceed in the analogous way. From 
equation (|40p . the axis data (j4"Tj) and the first condition from (fit)]) , we first arrive at the 
equation 

R(z)R*(z) [R{z)R * { ~ z) + R *W R (- Z ^ = °' ( 45 ) 
which is satisfied by all a n identically. Then we get 

R*(a n )R(-a n ) . 

— — — + 1 = X(a n ) = —X(—a n ), 46 

R(a n )R*{-a n ) 

thus arriving at the required restrictions on X n that characterize the equatorially 
antisymmetric solutions; namely, 

X k = -X 2N . k+1 , k = l,...,N. (47) 

The above identification of the equatorially symmetric/antisymmetric vacuum 
soliton solutions, carried out directly in terms of the parameters entering the expression 
of the Ernst potential defined by formulae (I36al) - fl36fj) . permits the consideration of 
particular configurations with additional equatorial symmetries independently of the 
corresponding axis data, which sometimes turns out to be very advantageous during 
the treatment of specific problems. In the papers [H] [15], for instance, the N = 3 and 
N = 4 specializations of the relations fl39l) and fT44l) were successfully employed for the 
analysis of equilibrium states in the triple- and quadruple-Kerr solutions. 



5. Concluding remarks 

In this paper we have emphasized the equatorially symmetric and antisymmetric 
electrovac solutions constructed from rational axis data. Such solutions have direct 
relevance to multi-black-hole spacetimes and to the fields of astrophysical objects. 
When an electrovac solution does not arise from rational axis data, the fundamental 
theorem proved in our previous paper can be employed to establish whether or not the 
solution in question is equatorially symmetric/antisymmetric. For example, from the 
conditions ([1]) one can very easily establish the equatorial symmetry of the electrovac 
solution for a charged rotating Curzon mass p3]. We should also like to point out that 
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J. Sod-Hoffs, who is a student of V. S. Manko, is currently working on the derivation of 
all the metrical functions for two of our antisymmetric solutions. He intends to submit 
for publication a paper concerning his analysis of the equilibrium conditions as well as 
the physical interpretation of these solutions. 
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